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We discuss processes leading to two photon final states in peripheral heavy ion collisions at RHIC. 
Due to the large photon luminosity we show that the continuum subprocess 77 -^ 77 can be observed 
with a large number of events. We study this reaction when it is intermediated by a resonance made 
of quarks or gluons and discuss its interplay with the continuum process, verifying that in several 
cases the resonant process ovewhelms the continuum one. It is also investigated the possibility of 
observing a scalar resonance (the a meson) in this process. Assuming for the a the mass and total 
decay width values recently reported by the E791 Collaboration we show that RHIC may detect 
this particle in its two photon decay mode if its partial photonic decay width is of the order of the 
ones discussed in the literature. 
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I. INTRODUCTION 

In the Rclativistic Heavy-Ion Collider (RHIC) operat- 
ing at the Brookhaven Laboratory beams of heavy-ions 
are coUiding with the main interest in the search of a 
quark-gluon plasma in central nuclear reactions. In ad- 
dition to this important feature of heavy- ion colliders, 
peripheral collisions may give rise to a huge luminosity of 
photons opening the possibilities of studying two-photon 
and other interactions as discussed by several authors pi . 

In order to study peripheral 7 — 7 heavy-ion collisions 
we have to pay attention to two important facts, which 
are fundamental to determine the kind of physics that 
is truly accessible in these interactions. The first one is 
to avoid events where hadronic particle production over- 
shadows the 7 — 7 interaction, i.e., events where the nu- 
clei physically collide (with impact parameter b < 2Ra, 
Ra being the nuclear radius) are excluded from calcu- 
lations of the usable luminosity. This exclusion has been 
discussed by Baur O] and can be handled as described 
by Cahn and Jackson H . The second point to remember 
is that the photons will carry only a small fraction of the 
ion momentum, favoring low mass final states. Therefore, 
hadronic resonances with masses up to a few GeV can be 
produced at large rates [Q,||, as well as any two-photon 
process not leading to extremely heavy final states. 

One of the purposes of our work was motivated by 
the possibility of measurement of the continuum process 
77 ^ 77 , see Fig.(||), as suggested in Ref. Q, as well as 
the resonant one 77 — *■ i? — > 77 which was predicted in 
Rcf . H as a clear signal for resonances made of quarks or 
gluons. The continuum 77 — > 77 reaction is interesting 
per se but is also a background for the resonant process. 
This last one is quite important to be observed because 
it involves only the electromagnetic coupling of the reso- 



nance. Its knowledge with high precision is very useful, 
for instance, to unravel the possible amount of mixing 
in some glueball candidates 0, complementing the in- 
formation obtained through the observation of hadronic 
decays. Another interesting study is the possible produc- 
tion of a light scalar meson (a) . 




FIG. 1. Diagram for 77 fusion in a peripheral heavy-ion 
collision. The blob represents a continuum or resonant process 
leading to a two-photon final state. 

The E-791 experiment recently reported that the a me- 
son has been observed in its tttt decay mode H, which 
is the same final state questioned over many years as a 
possible (or not) signal for this broad resonance. We dis- 
cuss the production of this elusive meson in peripheral 
heavy-ion collisions. However, we focus attention on the 
observation of its 77 decay. We believe that the observa- 
tion of this decay could settle definitively the existence 
of this particle as well as of its quark (or giuon) content. 
Although its branching ratio into photons is quite small 



compared to the main one (tttt), we will argue that the 
high photon luminosity in peripheral heavy-ion collisions 
may allow its observation, as long as it has the very same 
properties predicted in Ref. |Q. 

Photons with invariant mass of 0(1 GeV) are also one 
possible signal for the quark-gluon plasma, therefore, the 
detectors at RHIC have high efficiency to detect these 
photons, which is a welcomed capability in regard to the 
final states that we discuss in this work. The background 
for the two-photon initiated process were discussed by 
Nystrand and Klein M . Double-Pomeron exchange may 
be a background source Q. However, they are negli- 
gible after imposing the peripheral condition b > 2Ra 
as long as we operate with very heavy ions ^^. The 
cuts to remove the background are discussed in Ref. [^, 
i.e., the final state must have: a) low multiplicity, b) 
small summed transverse momentum and c) be centered 
around rapidity with a fairly narrow width of center- 
of-mass rapidity {i/cm ^ 1 ~ 2). These cuts are easily 
implemented when we analyze the reaction proposed in 
this work: ZZ —^ ZZjj. 

The remainder of the paper is organized as follows. In 
the second section we introduce the photon distribution 
in the ion that will be used in our calculations. Section 
III contains the discussion of the 77 — > 77 process and 
the cuts that we will adopt for this reaction. In section IV 
we present values for the production of resonances made 
of quarks (or giuons), and section V contains a discussion 
about the possibility of observing the a meson in its two 
photon decay. The last section contains our conclusions. 



II. PHOTON DISTRIBUTION IN THE ION 

The photon distribution in the nucleus can be de- 
scribed using the equivalent-photon or Weizsacker- 
Williams approximation in the impact parameter space. 
Denoting by F(x)dx the number of photons carrying a 
fraction between x and x -I- da; of the total momentum 
of a nucleus of charge Ze, we can define the two-photon 
luminosity through 
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dx 



F{x)F{tIx), 



(1) 



where r = s/ s, s is the square of the center of mass 
(c.m.s.) system energy of the two photons and s of the 
ion-ion system. The total cross section of the process 

ZZ — !■ ZZ^^ is 



cr(s) = / dT—a{s), 



(2) 



where a{s) is the cross-section of the subprocess 77 — > 
77. 

There remains only to determine F(x). In the liter- 
ature there are several approaches for doing so, and we 



choose the conservative and more realistic photon distri- 
bution of Ref. 1^ . Cahn and Jackson [gj , using a prescrip- 
tion proposed by Baur [Q| , obtained a photon distribution 
which is not factorizable. However, they were able to give 
a fit for the differential luminosity which is quite useful 
in practical calculations: 
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where z = 2MRy/T, M is the nucleus mass, R its radius 
and ^(z) is given by 



^(z) = ^A,e-''■^ 



(4) 



which is a fit resulting from the numerical integration 
of the photon distribution, accurate to 2% or better for 
0.05 < z < 5.0, and where Ai = 1.909, A2 = 12.35, 
A3 = 46.28, 61 = 2.566, 62 = 4.948, and ^3 = 15.21. For 
z < 0.05 we use the expression (see Ref. S) 
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(5) 



The condition for realistic peripheral collisions {bmm > 
Ri + R2) is present in the photon distributions showed 
above. 

The two photon production in heavy ions collisions 
can also be intermediated by the diffractive subprocess 
PP — > 77, where P is the Pomeron. But, as was shown 
by some of us in Ref. pOl, when the cut in the impact 
parameter is introduced the double-Pomcrom exchange 
is not important for very heavy nuclei. 



III. THE CONTINUUM REACTION 7 + 7^7 + 7 

The reaction 77 — *■ 77 (given by the fermion box dia- 
gram) was first calculated exactly by Karplus and Neu- 
man using the usual Feynman techniques pTJ ] . Some time 
later this process was calculated by De Tollis using the 
Mandelstam representation and dispersion relations p2| . 
They encountered relatively simple helicity amplitudes 
for this process. Of the sixteen helicity amplitudes, due 
to symmetry properties, the number of independent am- 
plitudes will be only five 0,II3], that may be chosen 

to be M++++, M++ , M+-+-, M+ + and M+++_. 

Where the + or — denotes the circular polarization val- 
ues -fl and —1. The remaining helicity amplitudes may 
be obtained from parity and permutation symmetry. 

Of these five helicity amplitudes, three are related by 
crossing, hence it is sufficient to give just three, that in 
agreement with Ref. |13,M are 
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where r^, St and t( are related with the standard Man- 
delstam variables s, t, and u by 
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ruf is the mass of the fermion in the loop. 

The remaining thirteen amplitudes may be obtained 
from these by using the relations 

M+++- = M++-+ = M+-++ = M-+++ 
M+-+-{s,t,u) = M++++{u,t,s) 
M+ — +{s,t,u) = M++++(t,s,u) 

(7) 

where s, t and u are the Mandelstam variables. 

The transcendental functions B, T and / tha t ap pear 
in Eq.(^) were defined by Karplus and Neuman [II|, and 



B[x) 



1 arcsinh(\/— a;) — 1, 

X 



1 arcsin('\/a;) — 1; , 

X 



1 arccosh(Va;) - 1 \l , {^ < x) (8) 

X 2 \ X 



T{x) ~ [arcsinh(-\/— a;)]^, 
— — [arcsin(-ya;)]^, 



(a;<0) 
(0 < x < 1) 



— [arccosh(v^)]^ i7rarcosh(-ya;), (1 < x) (9) 



i{x.v) = i{y,x), 

j?e{/(a;, y)} = ^ReU ( ^±^\ + $ f ^±^ 
^ ^ ''^^-' 2a \ \a + b{x)) \a-b{x) 



$( °',\ 1-$ 



a + 6(a;) 



a-1 
a — 6(2:) 



+ {x ^y) 



Im{Iix,y)} = -—\ny{a + bix))\ (x > 1) 
Za 

^-^\nx{a + b{yyf, (y > 1), 
Za 



(10) 



where a = -y/l — (x + y)/xy, x, y — rt, st and it, and 



6(x) = ^/ 1 - ( - ) , when {x < 0) or {x > 1) 



when (0 < a; < 1) 



(11) 



The Spence function, <I>(z), is defined as ^{z) — 
/o lii(l ~ t)dt/t, and its properties can been found in 
Ref. II . 

The differential cross section of photon pair production 
from photon fusion, i.e. the box diagram, is 
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is the scattering angle, a is the fine-structure constant, 
qf is the charge of the fermion in the loop and the sum 
is over the leptons e and /i and the quarks u, d and s, 
which are the relevant particles in the mass scale that 
we shall discuss. Another possible contribution to this 
continuum process comes from pion loops, which, apart 
from possible double counting, were shown to be negli- 
gible compared to the above one [G^. The second sum 



is over the sixteen helicity amplitudes, Mx-^ 



'^2'^3^4 



where 



Ai and A2 correspond to polarizations of incoming pho- 
tons and A3 and A4 for the outgoing photons. The matrix 
elements sumed over the final polarizations and averaged 
over the initial polarizations is given by 
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We consider the scattering of light by light, that is, 
the reaction 77 — > 77, in Au-Au collisions for energies 
available at RHIC, ^/s = 200 GeV/nucleon. We checked 
our numerical code reproducing the many results of the 
literature for the box subprocess, including asymptotic 
expressions for the low and high energies compared to 



the fermion mass present in the loop, and t he p eak value 
of the cross section (see, for instance, Ref. p3[ ). 

In Fig. (H) the dependence of the ion- ion cross section 
with the cosine of the scattering angle 9, in the two pho- 
ton center-of-mass system, is shown for an invariant pho- 
ton pair mass equal to 500 MeV. It is possible to observe 
that the cross section is strongly peaked in the backward 
direction, but is relatively flat out to cos 9 w 0.4, where it 
starts rising very fast. It is symmetric with respect to 9 
and the same behavior is present in the forward direction. 

The photon pair production in ion collisions over an 
invariant mass range between 300 and 1000 MeV for dif- 
ferent angular cuts is shown in Fig.(||). It is possible 
to observe that the cross sections is strongly restricted 
when the angular cut is more drastic, the difference be- 
tween the cuts 10° < 6* < 170° and 60° < 6^ < 120° is 
quite large, the cross section is strongly suppressed for 
the second range. We will impose in all the calculations 
throughout this work a cut in the scattering angle equal 
to I cos 9\ = 0.5. This cut is conservative, but it will make 
possible to compare the cross section of the box diagram 
with rival processes, that will be discussed in the follow- 
ing sections, as well as it is enough to eliminate the effect 
of double bremsstrahlung (which dominates the region of 
I cos^l w 1). Finally, this kind of cut is totally consistent 
with the requirements proposed in Ref. P] 
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FIG. 3. Invariant mass distribution of photon pair pro- 
duction at RIHC via photon fusion through a fermion loop. 
Each curve represents a different cut in the scattering angle, 
6, the solid line is for 10° < < 170°, the dashed is for 
40° < 61 < 140°, and the dotted one for 60° < 6) < 120°. 



FIG. 2. Angular distribution of ZZ -^ ZZ-y-y scattering at 
an invariant mass of 500 MeV. The scattering angle 9 is in 
the photon pair center-of-mass system. 

The fermions that contribute in the box diagram are 
the leptons e and /x and the quarks u, d and s, which 
are important for the mass range that we are interested 
(heavier quarks will give insignificant contributions and 
the same is true for the charged weak bosons). We as- 
sumed for their masses the following values: m^ — 0.5109 
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FIG. 4. Invariant mass distribution of photon pair pro- 
duction at RHIC center of mass energy of y^ = 200 GeV for 
each fermion loop. The solid line denotes the differential cross 
section for all fermions, the dashed one is the electron loop 
contribution, the doted is the muon one, the dash-dotted is 
the u quark contribution. The result is obtained with the cut 
Icos^l < 0.5. 



The electron gives the major contribution to the to- 
tal result as can be seen in Fig.(Q). The second most 
important contribution is due to the muon, but it is at 
least one order of magnitude smaller than the electron 
one. The d and s quark contribution (up to 0(2 GcV)) 
are smaller due to their masses and charges, because the 
process is proportional to (g?)^ where g/ is their charge. 
Their contribution is also insignificant compared to the 
electron result. 

As discussed in Ref. Q the 77 scattering can indeed be 
measured in peripheral heavy ion collisions. The cut in 
the angular distribution gives back to back photons in the 
central rapidity region free of the background. However, 
as we shall see in the next section, there are gaps where 
the 77 — *■ 77 process is overwhelmed by the presence of 
resonances like the rj, rji and others. Even the broad a 
resonance could be of the order of the continuum pro- 
cess. Just to give one idea of the number of events, with 
a luminosity of 2.0 xlO^^ cm~^s~^ |8| and choosing a 
bin of energy of 200 MeV, centered at the energy of 700 
MeV (which is free of any strong resonance decaying into 
two-photons), we have 1532 events/year assuming 100% 
efficiency in the tagging of the ions and photon detection. 



IV. RESONANCES DECAYING INTO 
TWO-PHOTONS 

Photon pair production via the box diagram is a back- 
ground to 77 — > i? — > 77 process (or vice versa), both 
have the same initial and final states, and for this reason 
they can interfere one in another. Normally the inter- 
ference between a resonance and a continuum process is 
unimportant, because on resonance the two arc out of 
phase. 

The total cross section for the elementary subprocess 
77 — > i? — > 77 assuming a Breit-Wigner profile is 
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(13) 



'fl^ total 



the well pronounced peaks of the resonances 77 and 77'. 
We assumed for their masses the values of 547.3 MeV 
and 957.78 MeV, respectively, the 77 total decay width is 
equal to 1.18 KeV and the 77' one is equal to 0.203 MeV. 
Their partial decay width into photons are 0.46 KeV (7;) 
and 4.06 KeV (77'). 

In the same figure we can see the predicted cross sec- 
tion for the glueball candidates production in peripheral 
heavy ion collisions by double photon fusion. 
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where M is the energy of the photons created by the 
collision of the ions. T[R -^ 77) (= T^-y) and Ttotai are 
the partial and total decay width of the resonance with 
mass 77ifl in its rest frame. 

We are going to discuss only J = resonances made 
of quarks as well of gluons. The reaction 77 — > 7r° — > 77 
was already discussed many years ago [ [f7| , where it was 
claimed that the interference vanishes. This result was 
criticized by De ToUis and Violini |]1^, affirming (cor- 
rectly) that the interference exists. However, as we will 
discuss afterwards, off resonance the interference is negli- 
gible. If the interference is neglected, Eq.(|l^) can be used 
and we show in Fig.(ph the result for some resonance pro- 
duction (77, 77', 77(1440), /o(1710)), whose invariant mass 
of the produced photon pair is between 500 MeV and 
2000 MeV. For comparison we also show the curve of the 
continuum process. It is possible to see in that figure 
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FIG. 5. Invariant mass distribution of photon production 
(with the cut | cos 8\ < 0.5). The solid curve is for the box dia- 
gram, the dashed curves are due to the process 77 — » _R — > 77, 
where R are the pseudoscalars resonances 77 and rj' and the 
glueballs candidates 77(1410) and /o(1710). 

We restrict the analysis to the J = glueballs candi- 
dates 77(1440) and /o(1710). For the 7^(1440) we used 
the mass and total decay width values of Ref. [[l9[ , 
TOfl = 1405 MeV, Ttotai = 56 MeV, for the decay width 
into photons we use the value given in Ref. |g^, T^^ — 5.4 
keV. We see in Fig. (pi) that the peak for this resonance is 
of the same order of the continuum process. For the other 
glueball candidate, /o(1710), the peak is clearly above the 
background. For this one we assumed the values listed 
in Ref. ||l^ of mass and total width, 771^ — 1715 MeV 
Ttotai — 125 MeV, and for the two-photon decay width 
we adopted the value encountered by the ALEPH CoUab. 



F. 



21.25 keV. In all these cases the resonances 



can be easily studied in peripheral heavy ion collisions. 

Off resonance we can expect a negligible contribution 
for the process 77 — > i? — > 77 and consequently the same 
for its interference with the continuum process. How- 
ever, it is instructive to present a more detailed argu- 



merit about why the interference can be neglected. In 
order to do so we are obliged to introduce a model to 
calculate the amplitudes for the process 77 — > i? ^ 77. 
These amplitudes will be computed with the help of an 
effective Lagrangian for the pseudoscalar interaction with 
photons (the scalar case will be discussed in the next sec- 
tion), which is given by gpe\^j^ypF'^^^F"''<^p, where gp is 
the coupling of the photons to the pseudoscalar field $p, 
£\fj,up is the antisymmetric tensor and F^^ is the electro- 
magnetic field four-tensor. 

The amplitudes for 77 — > 77 intermediated by a pseu- 
doscalar hadronic resonance are ilql: 
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M++__ = -^{F(AsO + F{\tt) + F{Xrt)}, (14) 

where a is the fine-structure constant, A — {raf /miif' , 
and 



ruR \ mil 



(15) 



The presence of the fine-structure constant in Eq. ( |l4| ) 
is a consequence of the fact that the amplitudes M in 
these equations will be used in Eq.(12), so it is necessary 
to get the correct dependence of the partial cross section 
with this constant. 

A numerical evaluation of the cross section using 
Eq.([l4|) (for the same resonances present in Fig.(|^)) 
shows a totally negligible effect off resonance in compar- 
ison with the box contribution. On resonance the two 
processes are out of phase and the interference is ab- 
sent. We can now proceed with an argument showing 
that the interference is not important. Let us assume 
that off resonance the processes are in phase, and for a 
moment we forget the t and u channels contribution in 
Eq. (p^ . The s channel contribution can be written as 
M2/Xs — rn?jl + iTiirafi), and denoting the continuum con- 
tribution by Ml we can write the following interference 
term 



^2^,2 



-^ — ^[(ReMiReM2 + ImMiImM^) 
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(s - m\) + {ReMiImM2 - ImMiReM2)TRmR\. 

We can verify that the term proportional to s — to|j inte- 
grates to zero when integrated in a bin centered at to|j. 
With the second term the situation is different. If ImMi 
or ImM2 ^ (assuming ReMi and ReM2 7^ 0) then 
there is a nontrivial interference. However, since we are 
dealing with J = amplitudes, the only nonvanishing 



helicity amplitudes are those in which the initial helici- 
ties and the final helicities are equal. Inspection of the 

A/++++ and M+^ amplitudes of the box diagram (in 

the limit mf <^ Tap;) shows that they are purely real, and 
the same happens with M2 (obtained from the s channel 
contribution of Eq.(|l4|)), resulting in a vanishing inter- 
ference! Of course, this analysis is model dependent. In 
particular, at the quark level the coupling gp has to be 
substituted by a triangle diagram, which may have a real 
as well as an imaginary part (see, for instance, Ref. p3|). 
However, this coupling is real for heavy quarks and its 
imaginary part is quite suppressed if the resonance cou- 
ples mostly to light quarks {ruf <^ fUR, what happens for 
the resonances that we are considering in the case of the 
u and d quarks, for the s quark the suppresion is not so 
strong in the 77 case, but we still have an extra suppres- 
sion due to its electric charge). Finally, the interference 
does appear when we consider the amplitudes with the 
resonance exchanged in the t and u channels, but it is 
easy to see that they are kinematically suppressed and 
also proportional to the small value of the total decay 
width. If the total width is large (as we shall discuss in 
the next section) the interference cannot be neglected. 
Although the above argument has to rely on models for 
the low energy hadronic physics, we believe that the di- 
rect comparison between the resonant and the continuum 
processes, as presented in Fig.(@), is fairly representative 
of the actual result. 

In Table I we show the number of events above back- 
ground for the 77, rj' and /o(1710) which, as shown in 
Fig.(^) are clearly above the box contribution. In the 
case of the /o(1710), as well as in the a meson case to 
be discussed in the next section, the decay of the res- 
onance into a pair of neutral pions is present. A pair 
of neutral pions can also be produced in a continuous 
two photon fusion process. The rates for all the reac- 
tions discussed in this work can be modified if both neu- 
tral pions, no matter if they come from a resonance or a 
continuos process, are misidentified with photons. This 
accidental background can be easily isolated measuring 
its invariant mass distribution, and making a cut that 
discriminates a single photon coming from the processes 
that we are studying, from one that produced two neu- 
tral mesons (mostly pions) subsequently decaying into 
two photons. For example, in the case of the sigma me- 
son (see next section) each one of the neutral pions from 
its much large hadronic decay should be misidentified. 
These pions would produce pairs of photons with a large 
opening angle 0, where cos ((/'/2) = -^/l — Am'^/'m^. 
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42 



TABLE L Number of events/year above background for 
the r], r]' and /o(1710) resonances. 



However, the calorimeters already in use in many exper- 
iments are able to distinguish between these two and single 
photon events with high efficiency (see, for instance, Ref. fej])- 



V. CAN WE OBSERVE THE a MESON IN 
PERIPHERAL COLLISIONS? 

The possible existence of light scalar mesons (with masses 
less than about 1 GeV) has been a controversial subject for 
roughly forty years. There are two aspects: the extraction 
of the scalar properties from experiment and their underly- 
ing quark substructure. Because the J = channels may 
contain strong competing contributions, such resonances may 
not necessarily dominate their amplitudes and could be hard 
to "observe". In such an instance their verification would be 
linked to the model used to describe them. 

Part of the motivation to study the two-photon final states 
in peripheral heavy-ion collisions was exactly to verify if we 
can observe such scalar mesons in its 77 decay. Although this 
decay mode is quite rare, it has the advantage of not being 
contaminated by the strong interaction of the hadronic final 
states. In particular, it may allow to investigate the possible 
existence of the sigma meson. This meson is expected to 
have a mass between 400-1200 MeV and decay width between 
300-500 MeV, decaying predominantly into two pions. Of 
course, another decay channel is into two photons, with the 
background discussed in Section III. 

Recently the E791 Collaboration at Fermilab found a 
strong experimental evidence for a light and broad scalar res- 
onance, that is, the sigma, in the D'^ —^ -k'tv^tt'^ decay |7|. 
The resonant amplitudes present in this decay were analyzed 
using the relativistic Breit-Wigner function given by 

BW= - 



mg -|- imoT{Tn) 



with 



mass and total decay width, while we vary the partial width 
into two photons. 

For the a decay width into two photons we assume the 
values obtained by Pennington and Boglione, 3.8 ± 1.5 KeV 
and 4.7 ± 1.5 KeV ||l|, and the value of 10 ± 6 KeV |o). 
The results obtained are in Fig.(te|), the dashed curve is for 
F-y-, — 3.8 KeV, the dotted one is for r^.v — 4.7 KeV and 



the dotted-dashed one is for F-^ 



10 KeV. In all cases we 



used the mass and total decay width obtained by the E791 
Collaboration, m^ = 478 MeV and Ftotai = 324 MeV. In the 
same figure it can be seen the curve due to the box diagram. 
In all cases we used the angular cut —0.5 < cosS < 0.5. 

It has been verified in the case of tttt scattering that the 
use of a constant total width in the a resonance shape is not 
a good approximation [El|. In our case we will discuss the 
77 ^ 77 process above the two pions threshold where the 
peculiarities of the broad resonance, basically due to the a 
decay into pions, are not so important. Of course, another 
reason to stay above 300 MeV is that we are also far from 
the pion contribution to 77 scattering. In any case we also 
computed the cross section with a energy dependent total 
width F(m) ~ Fp (p*/pn) ' , which, as shown by Jackson 
many years ago tea], is more appropriate for a quite broad 
resonance. The net effect is a slight distortion of the cross 
section shape with a small increase of the total cross section. 
Since this one is a negligible effect compared to the one that 
we will present in the sequence we do not shall consider it 
again. 



where m is the invariant mass of the two photons forming a 
spin-J resonance. The functions '^ F are the Blatt-Weisskopf 
damping factors [p5|: "_F = 1 for spin particles, ^F — 
l/y/Y+Jrp*Y for spin 1 and ^F = 1/^/9 -I- 3(rp*)2 + {rp*Y 
for spin 2. The parameter r is the radius of the resonance (« 3 
fm) [g6| and p* — p* [m) the momentum of decay particles at 
mass m, measured in the resonance rest frame, pg — Po("^)) 
where ttiq is the resonance mass. The Dalitz-plot of the decay 
can hardly be fitted without a O^""" (a) resonance. The values 
of mass and total decay width found by the collaboration with 
this procedure are 478^23 ± 17 MeV and 324^^^ ± 21 MeV, 
respectively. 

We will discuss if this resonance can be found in peripheral 
heavy-ion collisions through the subprocess 77 -^ ct ^ 77. 
It is important to note that all the values related to the a, 
like mass or partial widths, that can be found in the literature 
are very different and model dependent. In particular, we find 
the result of the E791 experiment very compelling and among 
all the possibilities we will restrict ourselves to their range of 
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FIG. 6. Invariant mass distribution of photon production. 
The solid curve is the box diagram, the dashed curve is for the 
process 77 ^ cr — *■ 77 when the F-y.y = 10 KeV, the dotted 
one is for F^y^. — 4.7 KeV and the dotted-dashed one is for 
T^.y = 3.8 KeV. The angular cut is equal to —0.5 < cos 6 < 0.5 
in all curves. 



Up to now we assumed in this analysis only the Breit- 
Wigner profile for the a contribution and it is interesting to 
discuss the full process (including the t and u channels) within 
a given model to verify how good is our approximation. 

To compute the amplitudes for 77 -^ 77 including the 
presence of a scalar resonance we can assume the interaction 
QsF^" Ffiv'^s, where Qa is the coupling of the photons to the 
scalar field $s- These amplitudes are 
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+++ = -— ^F{\Tt) 



= 0, 



2-K 



F{Xst 



{F{\st) + F{Xtt) + F{\rt)}. 



(16) 



until a final understanding comes out, with the advantage 
that the final state is not strongly interacting. Note that in 
this modelling the rj meson will contribute to 77 -^ 77 in a 
small region of momentum, even so it has to be subtracted in 
order to extract the complete a signal. 



Note that this interaction violates unitarity and is appro- 
priate only at energies up to the a mass. 

In Fig.(M) we show the full cross section for the process 
77 —> 77 due to the box diagram and the existence of the 
scalar resonance a with amplitudes given by Eq.(|lq), com- 
puted with F-y-y = 4.7 KeV. Contrarily to the previous cases 
the cr is a quite large "resonance" and due to its large decay 
width off resonance the interference is not negligible (the ar- 
gument of the previous section obviouly fails). Moreover, the 
interference is destructive! In Fig.dTI) we can see the continu- 
ous curve which is due only to the box diagram, the dashed 
one is the cross section of a scalar resonance computed with 
a Breit-Wigner profile, Eq.([13), the dash-dotted curve is due 
only to the scalar contribution of Eq.(pii[), and the dotted one 
is obtained when we consider the total cross section for two 
photon production in peripheral heavy ion collisions, i.e., the 
box diagram amplitudes plus the helicity amplitudes listed in 
Eq.(|lM). The figure shows clearly that the above mentioned 
interference is destructive. 

There is an important point to comment about the above 
result. The effective Lagrangian model used to compute the 
a contribution to the photon pair production above M ~ 600 
MeV gives a larger contribution to the cross section than the 
Breit-Wigner one. We verified that the t and u channels pro- 
duce small variations in the result obtained when we consider 
only the s channel, but this one already grows as any nonuni- 
tary amplitude calculated in this approach, and it is unclear 
for us how far we can trust in this model. At small masses 
the cross sections also differ in their behavior with s. 

There are possible improvements that we may think of to- 
wards a more reliable computation of the a meson exchange 
amplitudes, but all of them are going to be plagued by uncer- 
tainties as we remarked in the beginning of the section. There- 
fore, we adopt the pessimistic view that the Breit-Wigner 
result is just a lower bound for our process, because effec- 
tive lagrangians tend to increase the contribution at large 
masses where the background vanishes rapidly, increasing the 
signal/background ratio. The Breit-Wigner profile does the 
opposite. The physical result is probably between both ap- 
proaches. From the experimental point of view we would say 
that the reaction 77 —> 77 has to be observed and any devi- 
ation from the continuum process must be carefully modeled 
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FIG. 7. Invariant mass distribution of photon pair produc- 
tion. The solid curve is due to box diagram only, the dashed 
one is due to the process 77 ^ cr ^ 77 in the Breit-Wigner 
approximation, the dash-dotted is the scalar contribution of 
Eq.(|lq), and the dotted one is due to the total process, com- 
puted with box amplitudes plus the scalar boson amplitudes 
of Eq.([Lq). In all cases Fj^, — 4.7 KeV, and the angular cut 
is equal to —0.5 < cos6 < 0.5 . 



We changed the values of the a mass and total width 
around the central ones reported by the E791 Collaboration. 
We do not observed large variations in our result, but noticed 
that it is quite sensitive to variations of the partial decay 
width into photons. It is interesting to look at the values 
of the significance which is written as La signail \J Cohack and 
characterizes the statistical deviation of the number of the 
observed events from the predicted background. The signif- 
icance as a function of the two photons decay width of the 
sigma meson, with mass equal to 478 MeV and total decay 
width of 324 MeV, is shown in Fig. (pi), were we used a lumi- 
nosity of £ = 2.0 X 10^® cm~^s~^ at RHIC and assumed one 
year of operation. The significance is above 2a 95% confi- 
dence level limit for two photon decay width greater than 4.7 
KeV while for a 'ha discovery criteria can be obtained with 
F^,^ > 7.5 KeV. If F^-y turns out to be smaller the numbers 
of Fig.(H) have to be scaled appropriately. 



VI. CONCLUSIONS 




FIG. 8. Significance as a function of decay width into 
two photons, T^-y, for a sigma meson with mass equal to 
478 MeV and total decay width of 324 MeV. The solid 
curve was obtained integrating the cross sections in the in- 
terval 438 < M < 519 MeV, the dashed one in the interval 
300 < M < 800 MeV. 



The numbers in Fig.M were computed with the signal 
given by the Breit-Wigner profile, and the background by 
the pure box diagram. The solid curve was obtained inte- 
grating the cross sections in the interval 438 < M < 519 
MeV, corresponding to the interval of mass uncertainty, while 
the dashed curve resulted from the integration in the interval 
300 < Af < 800 MeV. Note that there is no reason, a priori, 
to restrict the measurement to a small bin of energy. This 
choice will depend heavily on the experimental conditions. 
When the signal is extracted from Eq.(|l6|) and the cross sec- 
tions are integrated in the larger interval (300 — 800 MeV), 
we obtain values for the significance more than one order of 
magnitude above the curves of Fig.(H), while the significance 
is below the curves of Fig.(H) up to a factor 10 when the cross 
sections are integrated in the smaller interval of energy. The 
Breit-Wigner profile probably gives an average of the results 
that could be obtained within different model calculations. 
Therefore, for values of T^-, already quoted in the literature 
the sigma meson has a chance to be seen in its two photon 
decay mode. The discovery limits discussed above refer only 
to a statistical evaluation. Our work shows the importance of 
the complete simulation of the signal and background includ- 
ing an analysis of possible systematic errors that may decrease 
the significance. 



Peripheral collisions at relativistic heavy ion colliders pro- 
vide an arena for interesting studies of hadronic physics. One 
of the possibilities is the observation of light hadronic reso- 
nances, which will appear quite similarly to the two-photon 
hadronic physics at e^e~ machines with the advantage of a 
huge photon luminosity peaked at small energies |l|. Due to 
this large photon luminosity it will become possible to dis- 
cover resonances that couple weakly to the photons |5|. 

In this work we focus the attention at the peripheral re- 
action ZZ —> ZZyy. This process is important because it 
may allow for the first time the observation of the continuous 
subprocess 77 — > 77 in a complete collider physics enviro- 
ment. This possibility only arises due to enormous amount of 
photons carried by the ions at the RHIC energies. 

The continuous subprocess is described by a fermionic box 
diagram calculated many years ago. We computed the pe- 
ripheral heavy ion production of a pair of photons, verifying 
which are the most important contributions to the loop, which 
turned out to be the electron at the energies that we are work- 
ing, and stablished cuts that not only ensure that the process 
is peripheral as well as eliminate most of the background. Af- 
ter the cut is imposed we still have thousands of photons pairs 
assuming 100% efficiency in tagging the ions and detecting the 
photons. 

The continuos 77 — > 77 subprocess has an interesting inter- 
play with the one resulting from the exchange of a resonance. 
We discuss the resonance production and decay into a photon 
pair. This is a nice interaction to observe because it involves 
only the electromagnetic couplings of the resonance. There- 
fore, we may say that it is a clean signal of resonances made 
of quarks (or gluons) and its measurement is important be- 
cause it complements the information obtained through the 
observation of purely hadronic decays. It may also unravel 
the possible amount of mixing in some glueball candidates 
p. We discuss the interference between these process and 
compute the number of events for some specific cases. 

The possibility of observing resonances that couple weakly 
to the photons is exemplified with the a meson case. This me- 
son, whose existence has been for many years contradictory, 
gives a small signal in the reaction 77 -^ cr -^ 77. However, 
its effects may be seen after one year of data acquisition, pro- 
viding some clue about this elusive resonance. Using values 
of mass, total and partial widths currently assumed in the 
literature, we compute the full cross section within a specific 
model and discuss the significance of the events. Our work 
shows the importance of the complete simulation of the sig- 
nal and background of these processes including an analysis 
of possible systematic errors, indicating that two photon final 
states in peripheral collisions can be observed and may pro- 
vide a large amount of information about the electromagnetic 
coupling of hadrons. 
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